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ABSTRACT 

We compare the stellar wind torque calculated in a previous work (Paper II) to the spin-up and 
spin-down torques expected to arise from the magnetic interaction between a slowly rotating (~ 10% 
of breakup) pre-main-sequence star and its accretion disk. This analysis demonstrates that stellar 
winds can carry off orders of magnitude more angular momentum than can be transferred to the disk, 
provided that the mass outflow rates are greater than the solar wind. Thus, the equilibrium spin 
state is simply characterized by a balance between the angular momentum deposited by accretion and 
that extracted by a stellar wind. We derive a semi-analytic formula for predicting the equilibrium 
spin rate as a function only of the ratio of Af„/il/a and a dimensionless magnetization parameter, 
5* = B'^RKMi^Vcsc)^^ , where A/^ is the stellar wind mass outflow rate, Ma the accretion rate, B^, the 
stellar surface magnetic field strength, i?» the stellar radius, and Uesc the surface escape speed. For 
parameters typical of accreting pre-main-sequence stars, this explains spin rates of ~ 10% of breakup 
speed for M^/M^. ~ 0.1. Finally, the assumption that the stellar wind is driven by a fraction of the 
accretion power leads to an upper limit to the mass flow ratio of M„/Afa ^ 0.6. 
Subject headings: accretion, accretion disks — MHD — stars: magnetic fields — stars: pre-main- 
sequence — stars: rotation — stars: winds, outflows 



1. INTRODUCTION 

The slow rotation rates of low to intermediate mass 
(< SA^q) pre-main-sequence stars remains one of the 
most important aspects of star formation that has, so far, 
resisted a generally accepted explanation. By the time 
they beco me optically visible as T Tauri stars (TTSs; 
lJovlll945l ). approximately half of them are observed to 
rotate at approximately 10% of breakup speed (the "slow 
rotators"; e .g., Vogcl & Kuhi 1981; Bouvier et al. 1997; 
iRebull et al.ll2004HHerbst et al.ll2007D . This is a surprise 
because many TTSs (the Classical T Tauri stars; CTTSs) 
are actively accreting material from sur r ounding Kep- 
lerian disks (Lvnden-Bell & Prine:lc"l974'; 'Bertout et al.l 
[£988; Calvct & GuUbring 1998; MuzcroUc et al. 2001). 
At a typical accretion rate of A^fa ^ lO^^Af© yr~^, the 
angular momentum deposited by accreting disk mate- 
rial should spin up a CTTS to near breakup speed in ^ 
10^ y ears (jHartmann fc Stauffeil 119891; iMatt fc Pudritd 
|2007( ) . Since the accretion phase lasts for 10^ - 10 '' 
years ()Lvo fc Lawsonl 120051 : iJavawardhana et "all l2006f ) , 
since the stars accrete at much higher rates prior to 
the TTS phase, and s ince the stars are still contract- 
ing (|Rebull et al.ll200^ . an efficient angular momentum 
loss mechanism is required to explain the existence of the 
slow rotators. 

A few interesting and important ideas for explaining 
the TTS slow rotators have been developed over the 
last two decades. T hese have result e d in t he star-disk 
interaction mo del of i Ghosh fc LambI (Il978f), appli ed to 
CTTSs byESnigl ( 19911 ajid_see lCamenzindl ll990l). the 
X-wind model (Shu et al]|1994D . and the idea that stel- 
lar w i nds provide strong torques (iHartmann fc Stauffcr 



1989'; 'Tout fc Pringle' 199^; Taatz & Camenzmdi 11996 ; 



Ferrciraet al. 2000; Matt fc Pudritz 2005a). Although 
both have advanced our understanding of the magnetic 



star-disk interaction, neither the G hosh fc Lamb nor X- 
wind models are without problerns (iFerreira et al.l 120001 ; 
lUzdenskvl l2004 iMatt fc Pudriti |2005^, and~lhe idea 
that stellar winds are important has not yet been worked 
out in sufficient detail to compare to the other models. 

In [Matt & Pudritz (2005a, hereafter Paper I), we fur- 
ther explored powerful stellar winds as a solution to the 
angular momentum problem and suggested that a frac- 
tion of the accretion power provides the energy necessary 
to drive the wind. We showed that stellar winds are ca- 
pable of carrying off the accreted angular momentum, 
provided that M^/Ma. ~ 0.1, where Afw is the outflow 
rate of material that is magnetically connected to the 
star (the "stellar wind"). This analysis included a for- 
mulation for the stellar wind torque that contained the 
Alfven radius (rA), which is not easily determined a pri- 
ori in the wind, and the conclusions were based on a one- 
dimensional scaling estimate of this important physical 
quantity. Thus, while it is clear that accretion-powered 
stellar winds (APSWs) can in principle provide the nec- 
essary spin-down torque, this idea requires further devel- 
opment to produce a more detailed model . 

Toward this goal. lMatt fc Pudritd (|2008l . hereafter Pa- 
per II) used 2-dimensional (axisymmetric) magnetohy- 
drodynamic simulations to solve for ta and calculate re- 
alistic stellar wind torques for a range of parameters. In 
the present paper, we use the stellar wind solutions of Pa- 
per II to compare the stellar wind torque to the torques 
expected to arise from the star-disk interaction. Fur- 
thermore, we find new solutions for stellar spins, based 
upon torque balance between the accretion torque and 
the APSW spin-down torque. This paper begins with 
a brief description of the simulation results of Paper II 
(Sj2]). We then compare the stellar wind torque to the 
star-disk spin-down torque in section [3] and then to the 
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TABLE 1 
Fiducial Stellar Wind 
Parameters 



Parameter 


Value 


M* 


0.5 Mq 




2.0 Rq 


B, (dipolo) 


200 G 


/ 


0.1 


Mw 


1.9 X IO-^Mq yr-l 




0.222 


7 


1.05 



star-disk spin-up torque in scctionlH which contains spin- 
equihbrium solutions. Section [5] contains a summary and 
discussion. 

2. RESULTS OF STELLAR WIND SIMULATIONS 

This section contains a brief description of the simu- 
lation results of Paper II that we will use for our anal- 
ysis, and the reader will find details in that paper. The 
primary purpose of the simulations was to compute the 
spin-down torque on a star, due to the angular momen- 
tum outflow in a wind. We used numerical magnetohy- 
drodynamic simulations to directly calculate the torque 
Tw from steady-state, 2D (axisymmetric) winds from 
isolated stars. We adopted coronal (thermal-pressure- 
driven) winds as a proxy for the unknown wind driv- 
ing mechanism. In the simulations, the torque is en- 
tirely determined by the seven key parameters listed in 
table [T] These are the stellar mass, M,; stellar radius, 
i?, ; strength of the rotation-axis-aligned dipole magnetic 
field at the surface and equator of the star, i?,; spin rate 
expressed as a fraction of breakup speed, 

f = n,Rl'^{GNU)-^''', (1) 

where £7* is the angular spin rate of the star; mass out- 
flow rate in the stellar wind, Mw; ratio of the thermal 
sound speed to the escape speed, evaluated at the base 
of the wind (just above the stellar surface), Cs/Vesc] and 
adiabatic index, 7. 

Table [T] lists the value of each parameter adopted for 
a fiducial case. Paper II contained a parameter study in 
which each of the seven parameters were varied relative 
to the fiducial case, and 14 cases from the parameter 
study are listed in table . In each case, six of the 
parameters were held fixed at the fiducial value (as given 
in table [1]), and one parameter was varied as indicated 
by the first column of tabled] 

To compare with analytic theory, we also calculated 
the effective Alfven radius (^a), where the poloidal wind 
velocity equals the poloidal Alfven speed, using an ana- 
lytic formula for the stellar wind torque. 



(2) 



Since our simulations are multi-dimensional, we have 
used {r\), which is the mass-loss- weighted average of 

1/2 

r\. Hereafter, we'll refer to {rj^) generically as ta- 
Using the simulation result for r^, equation [2] defines the 

^ In this paper, we do not discuss the cases from Paper II that 
include a quadrupole magnetic field, nor the extremely slow rotator 
case (with / = 0.004). 



value of TA, which is tabulated for all cases in the second 
colmim of table O 

In this paper, we make use of the semi-analytic formu- 
lation for the Alfven radius from Paper II, 



MwVcs 



(3) 



where K and m are dimensionless constants fit to the 
simulation, and Vcsc — {2GM^/R^y^'^ is the escape speed 
from the stellar surface. Paper II showed that the values 
of if w 2.11 and m « 0.223 well-describe (to better than 
1%) the fiducial case and those eight other cases with 
variations on B^,, i?*, Af^, and Af*. Although this is only 
approximately valid for situations with different wind ac- 
celeration rates or different rotation rates (in which case 
the values of K and m are slightly different; see Paper II), 
the formulation of equation [3] serves well as an indication 
of the approximate dependence of the stellar wind on pa- 
rameters, which will be important for discussing a wide 
range of possible conditions. 

Th e form of equation jp is si milar to t hat derived by 
(e.g.) iPelletier fc PudritTijigQa and see iFerreiral 119971 ) 
for the general theory of centrifugally driven disk winds. 
The quantity in brackets measures the magnetization of 
the wind. By assuming that the Alfven speed Vr^A (at the 
Alfven radius) is directly proportional to Q^^ta, a relation 
of the kind give n by equation Q can be de rived (e.g., see 
equation 2.27 of IPelletier &: Pudritdll99^ . In that case, 
the value of the index is to = 1/3. While this value is not 
far from the results of our numerical simulations, the dif- 
ference is significant. One key reason for this may be that 
disk winds are in the regime of so-called fast magnetic ro- 
tators, whereas the rather slowly rotating TTS are either 
slow magnetic rotators (where wind-driving forces dom- 
inate over centrifugal ones) or are intermediate between 
these two regimes (see Paper II). 

For the discussion that follows, it is useful to highlight 
how the lever arm (va) and wind torque responds to 
changing the mass load (the mass loss rate) of the wind. 
The fact that the Alfven lever arm in a hydromagnetic 
wind gets smaller as the mass load of the outflow in- 
creases, as is seen in equation ([3]), seems to suggest that 
the wind would become ineffective. This is certainly not 
true however, because equation ^ assures that an in- 
crease in wind mass loss rate leads to a net increase in 
the torque that the wind exerts upon the star (the net 
wind torque scales as A^T^"^™). This is the basic reason 
why, by having an outflow rate that is a substantial frac- 
tion of the accretion rate, an accretion-powered stellar 
wind can be effective in countering the accretion torque. 

It is our goal here to compare the stellar wind torque 
to the torque expected to arise from the star-disk interac- 
tion, and the latter has only been determined thus far for 
a dipolar geometry. So we only consider here the cases 
from Paper II with a dipole magnetic field. We also adopt 
the following assumptions. Paper II indicated that the 
details of the wind driving have a relatively small, but 
not entirely negligible, effect on the stellar wind torque. 
In the absence of a detailed model for how APSWs are 
driven, we assume that the velocity profile of an APSW 
does not differ substantially from our simulations (Pa- 
per II), so that the calculated torques are valid. Secondly, 
Paper II considered winds from isolated stars. Here, we 
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TABLE 2 

Stellar Wind Alfven Radii and Comparison 
TO Star-Disk Spin-Down Torques 





rA/ -ft* 


fw/Tdsd 


Tw / fdsd 






I — fl 1 \ 

\P — u.ij 


(ft — n ni \ 


fiducial 


6.97 


59 


490 


/ = 0.2 


6.26 


23 


200 


/ = 0.05 


7.65 


140 


1200 


B, = 400 G 


9.55 


27 


230 


= 2 kG 


19.3 


4.6 


39 


low Mw" 


11.8 


17 


140 


very low Mw" 


16.7 


6.7 


57 


R» = 1.5 Rq 


5.96 


86 


730 


-R* = 3 i?0 


8.75 


34 


280 


M* = 0.25 Mq 


7.52 


49 


410 


M, = 1 M0 


6.42 


70 


590 


Cs/Vesc = 0.245 


6.64 


53 


440 


Cs/Desc = 0.192 


7.23 


63 


530 


7 = 1.10 


7.79 


72 


610 



" The mass outflow rate in the low and very low 
Mw cases is 1.9 X 10"^° and 3.8 X W'^^Mq yr'^, 
respectively. 

will use the computed torques to develop the APSW sce- 
nario, in which stellar winds are accompanied by disk 
winds, accretion flows, and the general star-disk inter- 
action (see figure 1 of Paper 1) . In reality, the accretion 
disk blocks a portion of the stellar wind, and it is not 
clear how much this will affect the stellar wind torque. 
For the present study, we will assume that the presence 
of the disk and accretion will not significantly influence 
the stellar wind torques as computed in Paper II. 

3. STELLAR WIND VS. STAR-DISK SPIN-DOWN TORQUE 

The magnetic interaction between the star and disk 
results in angular momentum transfer between the two. 
All models that calculate the torque on the star from this 
interaction are based on the framework constructed by 
iGhosh fc LambI ()1978l ). In this general model, some of 
the stellar magnetic dipole flux connects to the accretion 
disk and conveys torques between the star and disk. The 
net torque can be separated into a spin-up part that adds 
angular momentum to the star and a spin-down part that 
removes angular momentum from the star, giving it back 
to the disk. In the absence of a stellar wind, a spin-down 
torque only arises when there is a magnetic connection 
between the star and the region of the disk outside the 
corotation radius. 



(4) 

It is assumed that the disk is capable of transporting 
away the excess angular momentum it receives from the 
star. The goal of this section is to compare the stellar 
wind torque to the spin-down torque arising from the 
star-disk magnetic connection, to determine under which 
circumstances each of these torques may be important 
and aid in the angular momentum loss from the star. 

To c alculate the star -disk spin-down torque, Tdsd, we 
follow iMatt &: Pudritd (|2005b, hereafter MP05b), who 
formulated a Ghosh & Lamb type model that includes 
the effect of the opening of magnetic field lines via the 
differential rotation between the star and disk. In this 
case, Tdsd is calculated by considering only the magnetic 
flux that remains closed (connected), parametrized as 
having an azimuthal twist of less than a critical angle. 



ta n~^(7c). Here we adopt the value of 7c = I suggested 
by lUzdenskv et ahl (|2002l ). By combining equation (9) 
and (22) of MP05b with equation dJ), we find 

xW) 



Tdsd 



(5) 



where 

x(/3)^rMl + (l + 2(1 + /?)-^] (6) 
is a dimensionless function of the strength of the effective 
magnetic diffusion rate in the disk. This is parametrized 
by , which for a standard a-disk (jShakura fc SunvaevI 
Il973l ) is /3 = aft.(Pt'")"^, where a has its usual meaning, h 
is the disk scale height at radius r, and Pt is the turbulent 
magnetic Prandtl number^ . Small values of /3 correspond 
to strong couphng, and MP05b suggested that (3 ~ O.OI 
was appropriate for real disks. For strong coupling (small 
P), the magnetic field will be highly twisted azimuthally, 
leading to more open flux and a weaker star-disk spin- 
down torque. For small /3, xiP) ~ but xiP) ti^s a 
maximum value of 0.25 when [3=1 (e.g., see fig. 7 of 
MP05b). 

Equation ^ indicates that the star-disk spin-down 
torque is completely independent of the accretion rate. 
This torque only requires that there exists a Keplerian 
disk outside Rco, to which the star can connect, and it 
does not matter whether or not there is net accretion 
onto the star. The dependence on the stellar spin rate / 
is due to the fact that when the star spins faster, i?co is 
closer to the star, where the magnetic field is stronger. 

We compare the stellar wind torque computed in our 
simulations to the star-disk spin-down torque by list- 
ing the ratio T^/Tdsd in the last two columns of table 
m We consider both a case with j3 = 0.1, resulting 
in X ~ 0.0826, and a case with j3 = 0.01, resulting in 
X ~ 0.00980. It is apparent that the simulated stellar 
wind torques are tens to hundreds of times greater than 
Tdsd, the larger difference existing for smaller values of (3. 

Thus, for the simulated winds, we see empirically that 
the stellar wind is much more effective at spinning down 
the star than is the star-disk connection. This can be 
understood qualitatively as follows. In the stellar mag- 
netosphere, any torque on the star is primarily conveyed 
by the azimuthal twisting of its magnetic field. In the 
case of the magnetic field connecting the star to the 
disk, there is a limit to how much the magneti c field can 
be tw isted before the connection is lost (Uzd enskv et al.l 
I2OOI MP05b). In the case of a stellar wind flowing along 
the fleld, there is no such limit on the twist. The larger 
the mass outflow rate in the wind, the less capable is the 
magnetic field to keep wind material corotating with the 
star, and so the larger will be the twist of the field. 

To show the dominance of the stellar wind torque more 
generally and to identify the circumstances under which 
it might not dominate, we will use the semi-analytic for- 
mulation of equation ^ . By combining equations ([T]) — 
([3]), one obtains 



l/m-2 



2d3 



fBtR 



(7) 



for the stellar wind torque. At first it may seem unusual 
that T„ is weaker when the magnetic lever arm length, 

^ Note that this /3 has no relation to the usual "plasma beta" 
parameter that often appears in MHD studies. 
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Fig. 1. — The ratio of the stellar wind torque to the spin-down 
portion of the star-disk interaction torque versus the magnetic lever 
arm length in the stellar wind. When the ratio is much greater than 
one, the stellar wind is most important for angular momentum loss 
from the star. The lines correspond to equation l|8j, assuming a 
stellar spin rate of / = 0.1 and three different values of the factor 
x(/3) (see text), corresponding to /3 = 0.01, 0.1, and 1, as indicated. 
The values listed in table[2] obtained by comparing simulated wind 
torques to the analytic disk torques, are plotted as squares (for 
f) = 0.1) and triangles (/3 = 0.01). The figure indicates that, 
unless the magnetic lever arm length is very long (e.g., for very 
low stellar wind mass loss rate), the spin-down torque from the 
disk is negligible. 

TA, is larger (for fixed fB^Rl). However this simply 
indicates that the stellar wind torque increases with in- 
creasing Mw, as noted in section [2l Also, note that a 
weak /-dependence of ta is not characterized in equa- 
tion ([3]) (see Paper 11), so the exact dependence of the 
torque on / is not captured in equation ([7]). 

By combining equations ([5]) and ([7]), using K = 2.11 
and m — 0.223 (see S|2]), one finds 

^.6.oxio^f^y7^y7^V".(8) 

The lines in figure [1] show equation ^ for / = 0.1 and 
for three different values of (3. This includes a line for 
f3 = 1, which corresponds to the strongest possible star- 
disk spin-down torque (as discussed by MP05b). Even 
in this case, the fiducial stellar wind torque is ^ 20 
times stronger than the star-disk spin-down torque. For 
smaller, more realistic values of /3, the star-disk spin- 
down torque only becomes weaker, while the stellar wind 
torque is not affected. It is clear that for the parameters 
considered here, the angular momentum extracted by the 
stellar wind completely dominates over that which can be 
transferred from the star to the disk. 

The stellar wind torque becomes weaker relative to 
Tdsd when TA is larger (e.g., for smaller M^) or for 
more rapidly spinning stars (larger /). For the case of 
/3 — 0.01, favored by MP05b, Tw will be larger than Tdsd 
for a star with / = 0.1, as long as ta ^ 84i?». This 
limiting value is much longer than any of the lever arm 
lengths listed in tabled As an example, for all else be- 
ing equal to the fiducial case, equation ([3]) suggests that 



Tw > Tdsd, as long as A/w ^ 3 x yr^^. This 

limit is comparable to the solar wind mass loss rate. If 
the stellar dipole field strength is instead i3* = 2 kG, the 
limit becomes Mw > 3 x lO'^^M©. 

The squares and triangles in figure [1] represent the data 
from table [D Note that nine of the data points (for each 
P) match the line very well. This is expected since the 
results from these cases were used to obtain the value of 
K and m used in equation ([5]). There arc five cases (for 
each /3) that lie slightly off of the line. Three of them 
represent the last three cases listed in table [21 and the 
other two are the / = 0.2 and / = 0.05 cases (which, for 
the plot, we have scaled by a factor of //O.l, to take into 
account the spin dependence of equation ([8])). These five 
cases are not expected to match exactly since equation 
([3]) is not precise for cases with different wind driving or 
spin rates than the fiducial case (see Paper II). Thus, the 
scatter of these five cases around the line indicates a sort 
of uncertainty of the semi-analytic formula for the stellar 
wind, due to variations in the wind driving mechanism 
and stellar spin rate. It is evident from the figure that 
this uncertainty does not affect the main conclusion here 
that stellar wind torques dominate the spin down of the 
star. 

Thus, for the slow rotators (/ ~ 0.1), we conclude 
that a stellar wind will transport much more angular 
momentum from the star than will a magnetic connec- 
tion to the disk, as long as the stellar wind mass out- 
flow rate is substantially larger than the solar wind mass 
outflow rate. For the systems considered here, with 
Mw ~ 1O~^M0 yr~^, the stellar wind torque completely 
dominates over any other spin-down torque felt by the 
star. Since Tdsd is negligible, the only important torques 
on these stars are the spin-down torque from the stellar 
wind and the spin-up portion of the star-disk interaction 
torque. We compare these two torques in the following 
section. 

4. SPIN EQUILIBRIUM BY AN APSW 

Section [3] revealed that, for the slow rotators with sub- 
stantial stellar winds, the spin-down torque felt by the 
star from the star-disk interaction is negligible. Thus, the 
spin state of the star is characterized as a competition 
between the spin-up component of the star-disk interac- 
tion torque and the spin-down by the stellar wind. If 
a system's parameters are measured, the theory can be 
used to determine the net torque on the star. 

Given enough time (~ 10^-10^ yr for CTTSs), the 
stellar spin should approach an equilibrium spin state 
in which the net torq ue on the star is zero (se e, e.g., 
ICameron fc Campbelll [l99l lArmitage &: Clarkd [l9961. 
The v ariability observed in accreting systems (e.g., 
iHartmann 1997) suggests that spin equilibrium may only 
represent a time-averaged state, and the condition of 
net zero torque simply identifies where the net torque 
changes sign. In any case, it is instructive to examine 
the conditions of spin-equilibrium. 

In sections 14.11 - 14. 3[ we examine the expected spin- 
equilibrium state of the specific cases of stellar winds 
simulated in Paper 11. In section 14. 4( we use the semi- 
analytic formulation of equation [3] to make more general 
conclusions. 

4.1. Spin Equilibrium for Specific Cases 
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TABLE 3 



Spin-Equilibrium Results for 7c 


= 1 & /3 


= 0.1 


Case 


. 

A/w/Ma 


Rt/R* 


btate 




fiducicil 


0.43 


4.4 


2 


0.71 


f — 2 


0.21 


2.9 


2 


0.39 


/ = 0.05 


0.83 


5.9 


1 


1.2 


B* = 400 G 


0.23 


4.5 


2 


0.37 


= 2 kG 


0.057 


4.6 


2 


0.055 


low Mw 


0.15 


4.6 


2 


0.23 


very low Mw 


0.076 


4.6 


2 


0.094 


R* = 1.5 Rq 


0.58 


4.3 


2 


0.96 


R* = 3 Rq 


0.28 


4.5 


2 


0.46 


M, = 0.25 Mq 


0.37 


4.5 


2 


0.61 


M« = 1 M0 


0.51 


4.4 


2 


0.84 


Cs/Vesc = 0.245 


0.48 


4.4 


2 


0.87 


Cs/'Uesc = 0.192 


0.40 


4.4 


2 


0.61 


7 = 1.10 


0.34 


4.4 


2 


0.48 



TABLE 4 

Spin-Equilibrium Results for 7;. = 1 & /3 = 0.01 



Case 


M^/M^ 


Rt/R, 


State 


^00 


fiducial 


0.44 


4.6 


2 


0.73 


/ = 0.2 


0.21 


2.9 


2 


0.39 


/ = 0.05 


0.83 


5.9 


1 


1.2 


B* = 400 G 


0.24 


4.6 


2 


0.39 


= 2 kG 


0.057 


4.6 


2 


0.055 


low M„ 


0.15 


4.6 


2 


0.23 


very low Mw 


0.077 


4.6 


2 


0.096 


iJ* = 1.5 Rq 


0.59 


4.4 


1 


0.98 


R* = 3 Rq 


0.28 


4.6 


2 


0.46 


M, = 0.25 Mq 


0.38 


4.6 


2 


0.63 


M, = 1 Mq 


0.52 


4.6 


2 


0.86 


Cb/vbbc = 0.245 


0.49 


4.6 


2 


0.88 


Cs/Vesc = 0.192 


0.41 


4.6 


2 


0.62 


7 = 1.10 


0.35 


4.6 


2 


0.49 



The spin-up portion of the star-disk interaction torque 
comes primarily from the accretion of material from the 
innermost part of the disk onto the star. In this section, 
we examine some specific cases of spin equilibrium, by de- 
termining under what conditions this spin-up torque bal- 
ances the stellar wind torques for the simulations listed 
in tabled 

When a star's magnetic field is strong enough, it will 
disrupt the Keplerian disk at some radius, i?t, the disk 
truncation radius. We calculate the location of i?t using 
the method and equations contained in the Appendix, 
which follows MP05b. From i?t, accreting material is 
channeled by the magnetic field to the surface of the star. 
There is a torque associated with the truncation of the 
disk and the accretion of material from i?f This torque, 
hereafter the "accretion torque," is given by MP05b as 



Ik 



1/2 



(9) 



where Ma is the mass accretion rate onto the stellar sur- 
face and k is the norm alized radius of gyration of the star 
(we assume ~ 0.2: lArmitage fc "Clarke 1996). Equa- 
tion ^ assumes that all of the Keplerian specific angular 
momentum of the disk material near Rt is transferred 
to the star. This naturally follows frorn the dynami- 
cal truncation of the disk (e.g., lYil 119951 IWangl 119951 : 
MP05b) and is supported by numerical simulations (e.g., 
iRomanova et al.ll2002l : iLong et al1l2005D . It is clear that 
the accretion torque depends both on Ma and Rt ■ At the 
same time, as shown in the Appendix, the location of Rt 
itself depends on most of the parameters, including Ma 
and p. 

The spin equilibrium state is defined by 

Ta = -Tw. (10) 

Each of our wind simulation cases represents a specific 
set of values for r^, Mw, i?*, M*, i?,, and O*. For each 
simulation case, we used equation ([9|) and the method in 
the Appendix to determine the values of Rt and Ma such 
that the condition pO)) is satisfied. We consider both a 
case with /3 = 0.1 and a case with (3 = 0.01 (see 
The results, given as M^/M^ and i?t/i?,, are listed in 
the 2nd and 3rd columns of tables [3] (for (3 — 0.1) andS] 
(for /3 = 0.01). 



A comparison between tables [3] and 2] reveals that the 
disk magnetic coupling parameter [3 has little influence 
on the equilibrium values of Ma and Rt- This demon- 
strates that, although /? has a large influence on the 
(negligible) spin-down part of the star-disk interaction 
torque (as shown in Sj3]), [3 has very little influence on 
the spin- up part. 

For the specific cases of the simulated stellar winds, 
it is clear that the equilibrium spin state is character- 
ized by Mw/Ma of typically a few tens of percent. This 
ratio is smaller for cases with larger rp^ (e.g., for larger 
field strength or smaller M^). Thus, the cases listed in 
the tables confirm the general conclusion of Paper I, and 
represent valid torque solutions for the spin equilibrium 
state. 

4.2. Magnetic Connection State of the System 

As pointed out by MP05b (and see the Appendix), an 
accreting system may exist in a state where the stellar 
magnetic field connects to the disk outside the corotation 
radius, which they call "state 2," and which we implic- 
itly assumed in section [3l On the other hand, if the disk 
truncation radius is sufhciently smaller than Rco, the star 
can lose its magnetic connection to all but the very inner 
edge of the disk, which they call "state 1." The determi- 
nation of Rt is different in the two states. In the absence 
of a stellar wind torque, a star in spin equilibrium must 
be characterized by state 2 (MP05b). Thus, having Rt 
very clo se to i?ro is a requirement of the "disk locking" 
models (iKonigllTiggl lOstriker fc Shul[l995l : IWang 1995). 
By contrast, this is not a requirement of the APSW sce- 
nario. 

Following MP05b (using equation lAljl . we determined 
the magnetic connection state of the spin-equilibrium 
systems described above and listed this in the 4th col- 
umn of tables [3] and m Note that MP05b only consider a 
loss of magnetic connection via the differential twisting 
of field lines. The stellar wind should also influence th e 
connectedness between the star and disk ()Safiedfr998f l. 
but we do not attempt to quantify this here. 

Tables [3] and S] reveal that most (though not all) of 
the simulated cases are in a magnetic connection state 2, 
while in spin-equilibrium. A characteristic of this state is 
that Rt is very close to i?coj which is also evident in the 
tables {Rco/R* = 7.4, 4.6, and 2.9 for / = 0.05, 0.1, and 
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0.2, respectively). It appears that, when a stellar wind 
torque balances the accretion torque, it may be common 
(though not required) for the disc truncation radius to 
be close to the corotation radius (unless the spin rate is 
substantially less than / = 0.1). 

4.3. Accretion Power 

In the APSW scenario proposed in Paper I, the energy 
that powers the stellar wind ultimately comes from the 
energy released by the accretion process. In this section, 
we calculate what fraction of the accretion power would 
be required to drive the wind, in the specific cases for 
which we have determined the spin-equilibrium. 

In order to tabulate the accretion power, the precise de- 
tails of the complicated interaction between the star and 
disk are not important. The general behavior is that 
material from the Keplerian disk becomes attached to 
the stellar magnetosphere near Rt and eventually falls 
onto and becomes part of the star. Energetically, this 
can be treated as an inelastic process, wherein only the 
energy content before and after the interaction needs 
to be specified. Thus, the rate of potential energy re- 
lease is simply 1/2 M^v'^^^^l — R^/Rt). Note that as 
accreting material piles onto the stellar surface, there 
should be additional energy released as material either 
"sinks" into (convectively) or compresses the star. This 
is another potential energy source, but we neglect this 
here. The rate of (rotational) kinetic energy release is 1 /4 

MaWesc(^*/^t — k^P), where the last term assumes that 
accreting material eventually achieves the same specific 
angular momentum as the star. The difference in ther- 
mal energy density between material at the disk inner 
edge and material at stellar photospheric temperature is 
negligible compared to the potential and kinetic energy 
release. Thus, by summing the potential and kinetic en- 
ergies, the rate of energy release in the vicinity of the 
star is approximately 



lUstvugova et al.lll999t iKeppens fc Goedbloedll2000f) 



4 = \My,,, 



^ 2R, 2-' 



(11) 



Since there is an accretion torque on the star, some of 
this energy is added to the rotational energy of the star 
at a rate ri*Ta. The remaining energy (E'a — ri*Ta) is 
available to power other accretion-related activity'^. In 
particular, this remaining energy is responsible for pow- 
ering the observed excess continuum emission (s uch as 
the UV excess) and line emission (e.g., Konigi fl99l 
iCalvet fc Gullbrindfl998l: iMuzerolle et al.ll200li) . in ad- 
dition to driving an enhanced stellar wind (Paper I). 

Paper I proposed that a fraction e of this energy specif- 
ically powers the thermal energy in the stellar wind. Our 
simulated winds are thermally driven, but the wind driv- 
ing mechanism at work in real systems is still uncertain 
(see Paper II). Thus, we wish to calculate the power re- 
quired to drive the wind, in a generic form. For this, we 
simply calculate the total energy in the wind far from the 
star plus the potential energy required to lift the wind off 
the stellar surface. In this way, the power in a steady- 
state, 2.5D, MHD wind can be obtained by (see, e.g., 

^ We follow a similar derivation of the accretion power to that 
of Paper I, where equation (4) of that work corresponds to the 
remaining energy, i?a — f^tfa, and neglects terms proportional to 



E„ = 47ri?2 / pv^E'd{cos 9) 



-Mw«4c, (12) 



where 9 is the usual spherical coordinate and 



E' 



(13) 

2 47rp AnpVp 

In equation (fT^]) . we have neglected the thermal and 
gravitational potential energy, so the integral in equa- 
tion should be evaluated at large i?, where E' has 
reached an asymptotic value and these energies are neg- 
ligible. Thus, i^w represents the total power required to 
lift material off of the star, to accelerate it to the wind 
velocity, and to provide the magnetic energy content car- 
ried with the wind. 

For each of our simulated wind solutions, we evaluate 
the integral in equation at a radius of i? = 50i?*, 
where E' is within a few percent of its asymptotic value. 
The spin of the star does work on the wind at a rate 
rj^Tw. This represents the power injected in the wind 
by magnetocentrifugal processes. We find that the ra- 
tio n^Tvi/Eyj is 30% in the fiducial case. In most other 
cases, the value of this ratio falls the range 10-60%. This 
indicates that, as discuss ed by Paper I and II (and see 
IWashimi fc Shibata"l993). these winds are in a regime 
where the magnetocentrifugal effects are of nearly equal 
importance with the other source of wind driving. 

It is this other source of wind driving that we propose 
is powered by some fraction of the available accretion 
energy. We define this fraction as'' 



E^, — fi* 



(14) 



This represents the minimum fraction of the accretion 
power required to drive the stellar wind, since whatever 
mechani sm drives the wi nd will not itself likely be 100% 
efficient pecamplilll98l . 

In the last column of tables [3] and [H we list the value 
of for each case in spin-equilibrium. In one case 
(/ = 0.05), £00 is greater than 100%, indicating that 
there is not enough accretion power in the spin equilib- 
rium state to power the wind. This case is therefore not 
an acceptable solution for a system in spin-equilibrium 
by an APSW. All of the other cases have eoo < 1, and so 
they are energetically viable solutions. 

There is a relationship between Eoq, the observed ex- 
cess emission, and the inferred mass accretion rate. In 
particular, the accretion rates are typically determined 
by measuring excess emission and assuming that all of 
the a ccretion power is radiated fe.g. . ICalvet fc Gullbrind 
IT99I . In the APSW scenario, some of the accretion 
power drives the stellar wind, so only a fraction 1 — eoo 
of the accretion power can be radiated. This means 
that the true accretion rate (Ma) will be a factor of 
(1 — £00)""'^ larger than the observationally determined 
value. In this context, and since eoo is the miminum re- 
quired fraction to drive the wind, a value of eoo ^ 0.5 

* This fraction is a more general definition than e in Paper I, 
which assumes thermal wind driving. By contrast, eoo is the frac- 
tion of the accretion power required to explain the energy in the 
wind at large distances from the star, regardless of the driving 
mechanism. 
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(as for several cases listed in tables |3] and |4]) appears 
quite large. However, the observational determination 
of Ma is uncertain by a factor of several, as exempli- 
fied by the large rang e of m easurements compiled by 
iJohns-Krull fc GaffordI (|200^ . Thus, while it is clear 
that £00 < 1 is a hard upper limit, it is not yet clear how 
close to unity e^o can be. 

As expected, the cases with lower values of M^jM^ 
(i.e., cases with larger field strength or smaller Mw) re- 
quire a smaller fraction of the accretion power to drive 
the wind. The cases in the table suggest approximately 
that eoo ~ 1.6Mw/Ma. Thus, it appears that Af^/Ma < 
0.6 represents a hard upper limit for APSWs. 

It is important to note that the manner in which the 
accretion power transfers to the stellar wind is still un- 
specified in the APSW model. This will depend upon 
what is the wind driving mechanism, which is currently 
unknown. In reality, the physics of the energy coupling 
will likely determine the value of eoo, which effectively 
sets the value for M^jM^. Then, given enough time, the 
stellar spin rate will evolve toward the equilibrium value 
set primarily by M^/Ma and rp^jR^. Thus the spin- 
equilibrium state of the star is ultimately determined by 
the power coupling and magnetic properties, and more 
work is needed to take this further. 

4.4. Semi-Analytic Formulation for Spin Equilibrium 

In order to develop a more general, predictive theory, in 
this section we make use of the semi-analytic formulation 
of equation ([3]). As justified by the previous sections, we 
assume that the equilibrium spin rate of the star is simply 
determined by a balance between the spin-up torque from 
accretion and the spin-down torque from the stellar wind. 
Using equations H]) - Q, (HI), and (fTO)) . we can write the 
stellar wind equilibrium spin rate, expressed as a fraction 
of breakup spin, as 



/sw = K- 




-2m 



where 



* = 



mi 



(15) 



(16) 



is a dimensionless magnetization parameter^. Here, we 
have neglected the term proportional to k'^ f in equation 
([9|. since it is generally much smaller than the other 
term. Again, note that a weak /-dependence of ta is 
not included in equation ([3]), so the dependence of /sw 
on some of the parameters is not precisely captured in 
equation (see Paper II). 

Equation (|15[) includes a dependence on the truncation 
radius of the disk, i?f This location itself has a depen- 
dence on the other parameters, and the determination 
of i?t depends on the magnetic connection state of the 
system f ii4.21 MP05b). In general, Rt depends on but 
if Rt is close to the corotation radius, i?co, then Rt also 
depends on the stellar spin rate. We will consider two 
cases that are expected to bracket reality. 

^ The magnetization parameter ^ is related to ip used by MP05b 
(see their eq. 16) by a constant factor, ip = 2'^/^^. 



The first case is one in which the system is in state 
1 as defined by MP05b. Here, Rt does not depend 
on the stellar spin rate, and i t is simply proport ional 
to the original calcu l ations by iLamb et al.l (|l973f) and 
[Davidson fc Ostrikerl (|1973[ ) used in most Ghosh & Lamb 
type models. Thus, using equation (jA2[) for Rt, adopting 
7c = 1, and plugging in to equation psp . one finds 




^l/7-2m_ 



(17) 



This is the predicted equilibrium spin rate when the trun- 
cation radius is significantly smaller than the corotation 
radius (i.e., in magnetic connection state 1). 

The second case to consider is where Rt « i?coi which 
is the requiremen t of all disk-locking models ()Shu et al.l 
119941: lWanslll995l MP05b). Using equation (U) and set- 
ting Rt = i?co in equation (|15p . one finds 



(3-6m)/4 



(18) 



This is the predicted equilibrium spin rate when the disk 
truncation occurs very close to i?co. 

Which case is more appropriate? The first case is ex- 
pected to occurs for relatively small values of ^E" and low 
spin rates / (MP05b), and the opposite is true for the 
second case. For the parameter space we have consid- 
ered thus far in this work, we found in section 14.21 that 
most (though not all) of the cases are expected to have 
Rt/ Rco near unity. Thus, while it is not a formal require- 
ment of APSW, it may often be the case that Rt « Rco 
for systems in spin equilibrium, and we will focus on this 
second case for the remainder of this work. 

Figure [2] shows the predicted spin rate of this case 
(eq. [T5)l versus for many different values of the ra- 
tio Mw/Afa. The results of section l4?3l indicate that the 
accretion power is only capable of powering a wind with 
Afw/Ma < 0.6, and this is a hard upper limit. This "for- 
bidden" region of the /-^ space is indicated in figure 

m 

The plot also shows the simulation results (squares). 
The value of "if for each case is determined mostly by 
input parameters but also by Ma. The latter was set (in 
section 14.11) by the condition that the equilibrium spin 
rate was equal to the value of / used as the simulation 
input parameter. The squares indicate the range over 
which equation is shown to be valid by the simula- 
tions. Also, the plot shows the specific case of = 2.11 
and m = 0.223. 

The results shown in figure[2]and equation (fT8|) provide 
the basis for predictions of the APSW model that can be 
observationally tested and constrained, and that can be 
compared to other models. Different theories predict dif- 
ferent power laws for / vs. Specifically, equation p8p 
predicts a power law index of « —0.33, the conditions 
appropriate for equation (fT7|) predict « —0.30, and an 
index of —3/7 is predicted by the disk locking models 
(e.g. , lKonig]|[l99ll : iShu et al.lfl99l . 



5. SUMMARY AND DISCUSSION 

In this work, we have further developed the accretion- 
powered stellar wind model proposed in Paper I, where 
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Fig. 2. — The equihbrium spin rate predicted by a balance 
between the spin down from a stellar wind and the spin up 
from accretion, versus the dimensionless magnetization parame- 
ter ^' = B^R1(M!^Vcsc)~^ ■ The solid lines show equation JTSl l 
for K = 2.11, m = 0.223 and several different possible values of 
M-w/Ma, as indicated. The squares show data from table [S] which 
indicates the range of parameters considered in our simulations 
and used to derive equation I ll8t . The shaded region corresponds 
approximately to where accretion power is not sufficient to drive 
the stellar wind. 

the stellar wind magnetic lever arm length, r\, was taken 
as a parameter to determine the stellar wind torque. We 
employed the simulation results of Paper 11 (see ^ to 
obtain stellar wind torques for several cases representa- 
tive of T Tauri systems. We examined the total torque 
on the star arising from the stellar wind plus the mag- 
netic interaction between the star and its accretion disk. 
Our results can be summarized as follows. 

1. We found that the spin-down torque from a stel- 
lar wind can be orders of magnitude stronger than 
the spin-down portion of the star-disk interaction 
torque, for slowly rotating stars with mass loss 
rates substantially larger than the solar wind out- 
flow rate (see @. This confirms the assumption of 
Paper I that the condition for net zero torque on 
the star (spin equilibrium) is simply determined by 
a balance between the stellar wind torque and the 
accretion torque. 

2. Using the computed stellar wind torques for several 
cases, we looked at the conditions for spin equilib- 
rium f ij4.ip . We found that a rotation rate of 10% 
of breakup speed typically requires M^/Ma. equal 
a few tens of percent, confirmin g the o riginal sug- 
gestion bv lHartmann fc Stauffeil ()1989[ ) that stellar 
winds may be capable of removing accreted angular 
momentum. 

3. For most cases in spin equilibrium, the disk trun- 
cation radius was very close to the corotation ra- 
dius, though this is not a general requirement of 
the APSW model (g^). 

4. Accretion power is generally sufficient to power a 
stellar wind that is capable of solving the angular 
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momentum problem ( !j4.3|) . as suggested in Paper 1. 
The energy requirements for most of the cases con- 
sidered here is relatively large, and more work is 
needed to further constrain the energy coupling. 

Under the assumption that the stellar wind is ac- 
cretion powered, the cases we examined suggested 
a hard upper limit of M^jM^ < 0.6. 



6. Finally, in section 14.41 we used the results from 
Paper 11 to derive a semi-analytic formulation for 
the equilibrium spin rate predicted by the APSW 
model. We found the that spin rate, expressed as a 
fraction of breakup speed, generally depends only 
on the two dimensionless parameters M^/Ma and 

The APSW model incorporates several previous ideas. 
As in all other models that emphasize the role of stellar 
magnetic fields, the interaction of the magnetized star 
with the disk leads to the truncation of the disk and ac- 
cretion of material along field lines onto the star. The 
APSW model adop ts the finding of the Ghosh fc Lamb- 
type models Ce.g.. [Ghosh fc LambI flOTSt IKonigll [l99ll : 
lArmitage fc Clarkg|1996D, which is also supported by nu- 
merical simulations (R omanova et all I2002t iLong et all 
2005), that the angular momentum of accreting material 
is transferred to the star. However, in contrast to the 
Ghosh fc Lamb-type models, we found that for slow ro- 
tators, any spin-down torque arising from the star-disk 
interaction is negligible (item 1 above). Instead, the stel- 
lar spin-up torque from accretion is counteracted by an 
accretion-driven stellar wind, which carries a comparable 
amount of angular momentum out of the system. 

Compared to all existing models, APSW is distinct 
in that it conceptually links the driving of the stel- 
lar wind to the energy released by the accretion pro- 
cess (via e). In other ways, the general picture of 
APSW is similar to other angular momentu m models 
that utilize winds. In partic ular, the X-wind ()Shu et al.l 
11994"; 'Ostriker fc Shul ll995D . the Reconnection X-wind 
^jgrrcir a et al. 200(j)7and other wo r ks considering stellar 
wind s (jHartmann fc Stauffeil 119891 : iPaatz fc Camenzindl 
Il996l ) all find that, in order to carry away significant an- 
gular momentum, the mass outfiow rate needs to be of 
the order of 10% of the accretion rate (item 2 above). 
Except for the X-wind, in all of the above scenarios, the 
outflow is magnetically connected to the star, and thus 
extracts angular momentum directly from the star. By 
contrast, the X-wind outflow is magnetically connected 
to the disk. Furthermore, the X-wind is unique in that it 
assumes that the accretion of material does not deposit 
angular momentum onto the star. 

As mentioned in item 3 above, for most of the specific 
cases of our simulated winds, we found that in the spin 
equilibrium state, the truncation radius was very close to 
the corotation radius. This is similar to the prediction of 
the disk locking models, which includes both the Ghosh 
fc Lamb-type models and the X-wind. However, in con- 
trast with the disk locking models, it is not a requirement 
of APSW that i?t be close to R^o- Measure ments of the 
location of the inner edge of the gas disk (jNaiita et al.l 
f200l [Canl[2007f) suggest that i?t/i?co is typically ~ 70%. 
We leave a more detailed comparison between models for 
future work. 
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There is observational evidence that the outflow rates 
from accreting young stellar systems are of the or- 
der of 10% of t he accretion rates (e.g., iHartigan et al.1 
1997„) and are therefore accretion powered 
II99OI ). However, it appears that a large 



fl995t ICalvetl 



(Cabrit et al 



fraction of this flow (which is usually probed by for- 
bidden emission coming from large spati al scales) orig- 
inates in the disk, rather than the star (|Ferreira et al.l 
[2()06). It is not yet clear how much of the total ob- 
served flow may originate in a stellar wind. There is 
some evidence specifically f or stellar winds from CTTSs 
feeristain et al. 2001; Edw ards et al.ll2003l: [Pupree et all 
[2005i , Edwards et al..,2006 : iKwan et al.ll2007| ). as distinct 
from disk winds, and that these are accretion powered 
(e.g.. lEdwards et al.|[2003l l2006'). b ut the mass outflow 
rates are not yet well constrained (Dupree et al.l [20051 ). 
Additional work constraining the value of M^jM^^ the 
stellar wind driving mechanism, and the stellar magnetic 
field strength and geometry will help to provide stringent 
and quantitative tests for the APSW model. 

The predictions of the spin equilibrium state can 
also be checked observationally (|Johns-Krull fc GaffordI 
I2OO2I ). This will likely require large samples of stars, due 
to large uncertainties in measured para meters, and since 
intrinsic variability in real systems ( HartmannI Il997l 
e.g.,) may only allow a spin equilibrium state to be 
achieved in a time-averaged sense. The Ghosh & Lamb, 
X-wind, and APSW models all predict an equilibrium 
spin rate that depends on but of these three, only the 
APSW model contains an additional dependence on the 
stellar wind mass outflow rate. For the power law fits to 
the simulations of Paper II, the APSW spin equilibrium 
predicts a slightly different power-law of spin vs. ^' than 
the other models ( §4.4p — though the exact dependence 
of the stellar wind torque has not been determined for 



all parameters. 

In this series of papers, we have focused on the global 
problem of calculating the magnitude of stellar wind 
torques and comparing them with other torques acting on 
accreting stars. In order to refine the APSW model fur- 
ther and make the predictions more precise, more work 
is required. In particular, it is not yet clear how the pres- 
ence of an accretion disk will influence the stellar wind 
torque, and conversely, how a stellar wind may influence 
the accretion process. Also, although it is clear that there 
is enough accretion energy to power the stellar wind, it 
is still not known what actually drives the stellar wind 
and how the accretion power may transfer to it. We sus- 
pect that a strong flux of hydromagnetic waves can be 
excited near the base of the accretion shock and can tap 
the energy released there, which may provide an efficient 
driver for the APSW. We defer a rigourous investigation 
of the APSW driving mechanism to future work. 
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APPENDIX 

DETERMINATION OF THE DISK TRUNCATION RADIUS 

We follow MP05b to calculate the location of the disk truncation radius, i?t, and the reader will find details in that 
paper. For convenience, we list the relevant equations here. As in section [Sj we adopt 7c = 1. 
To determine i?t, we first find the magnetic connectivity state using the criterion 

/ < (l-/3)(23/2*)-3/7^ (Al) 

where 4' is defined by equation p6|l . If condition (|Aip is satisfied, the system is in "state 1." Otherwise, it is in "state 
2." In state 1, we determine the truncation radius using 

i?t = (2^/2^)^/'^^*. 

In state 2, we determine the truncation radius by solving 



(A2) 



^co 



-7/2 



1 



Rco 



3/2' 



23/2^-- 



-7/3 



(A3) 
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